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We compare observables to the finite size scaling hypothesis in SU(2) lattice gauge theory with 
two Dirac fermions in the adjoint representation. The fits that we obtain yield an estimate of the 
anomalous mass dimension that is consistent with four loop perturbation theory: 7 = 0.51 ± 0.16, 
with the error due to systematic uncertainties in the finite size scaling analysis. The result is 
somewhat larger than one Schrodinger functional study (by 1.3a) but consistent with another. 
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I. INTRODUCTION 

In technicolor models, the Higgs mechanism occurs 
through condensation of new fermions that are subject 
to a gauge interaction that is strong at the TeV scale 
P, Walking technicolor is a version of this theory 
that can suppress flavor-changing neutral currents by 
raising the extended technicolor scale, while still having 
phenomenologically acceptable Standard Model fermion 
masses, due to condensate enhancement [l-Q- Higher 
representations of the gauge group are believed to avoid 
problems with the S-parameter, i.e. electroweak precision 
constraints [13, El • All of this has motivated the study 
of Minimal Walking Technicolor (MWTC) [13], which is 
SU(2) gauge theory with two Dirac fermions in the ad- 
joint (triplet) representation. 

In order to study technicolor proposals nonperturba- 
tively and from first principles, several groups have been 
using the techniques of lattice gauge theory; see the re- 
view [131 ] and references therein. One of the key ques- 
tions is whether the theory "walks" (very slow running 
of the coupling) or is attracted to an infrared fixed point 
(IRFP). An important quantity that can be computed in 
the process of answering this question is the anomalous 
mass dimension 7, which needs to satify 7 ~ 1 in order for 
the standard walking technicolor picture to succeed. (Al- 
ternatives such as "ideal walking" are now being inves- 
tigated as improvements over the standard picture [Til].) 
One of the ways in which the lattice community has com- 
puted 7 is through the Schrodinger functional method. It 
was employed for SU(3) gauge group with sextet fermions 
in [H| and for MWTC in 0,ll3- Analysis of the distri- 
bution o f eig envalues of the Dirac operator has also been 
used PHI}. 

An alternative approach is to compare observables 
computed in lattice gauge theory (e.g., meson masses, 
the "pion" decay constant) to the finite size scaling (FSS) 
hypothesis. If the theory is indeed driven to an IRFP, 
then the data on observables should fit the FSS hypothe- 
sis. Previous studies of FSS in lattice technicolor include 
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(19l |21|-|23j. Fits to the conformal hypothesis that as- 
sume a specific form of the FSS function include 0, HH , 
where the infinite volume hyperscaling relation is im- 
posed. More general forms of the FSS function have also 
been considered recently by the authors of |24|, with the 
result that for these forms the resulting conformal hy- 
pothesis for SU(3) gauge group and 12 fundamental fla- 
vors has a low de gree of confidence in fitting the data 
[26j|. By contrast, |23l ] advocates an approach that does 
not impose a specific form on the FSS function; this is 
one of the FSS methods used in the earlier work (l9j . In 
this letter, we apply this method to the case of MWTC 
in order to extract an estimate of 7 under the assumption 
that an IRFP exists. 



II. FINITE SIZE SCALING 

In the scaling regime, the correlation length will have 
an asymptotic behavior dependent on the fermion mass 
m with exponent y m : 



-l/Vn 



(2.1) 



This exponent is related to the anomalous mass dimen- 
sion evaluated at the IRFP: 



Vm = 1 + l(9*) 



(2.2) 



It is a general consequence of the renormalization group 
equations that the correlation length £l at finite size L 
is given by a scaling function of the infinite volume cor- 
relation length £ relative to L: 



i L jL = F{£/L) 



(2.3) 



Thus we obtain the FSS formula in terms of fermion 
mass: 



&/L = /(mi"") 



(2.4) 



Corrections to scaling will be an important consideration 
for us. This translates into a correction that is apprecia- 
ble for small L, with an exponent u>: 



£l/L = /(mP») + L-"g(mLy™) 



(2.5) 
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This form has also been considered in [26[; there it was 
pointed out that fitting data to such a hypothesis would 
require an extensive and highly accurate study. For us 
the main use of this equation is just that the scaling 
violations are largest for the smallest values of L. We 
use this as an interpretation of data on small lattices 
that does not fall on a scaling curve. Our present study 
is not extensive enough to fit to this more general form 
and extract uj. Below, we will consider £x = 1/M or 
I/Ztt, where M is a meson mass. 



III. FITTING METHOD 



The method described here seeks to optimize y m such 
that all the data falls on a scaling curve. It is due to [I?} 
and was used in [l9|, HH . For each L we have a data set p. 
We use this to obtain a fit f p . The types of fit functions 
that we consider will be described below. We then use 
this fit function on the other values of L, which we label 
as L .j. 



We minimize the following function with respect to y n 



v^- EE E - m^)) 2 (3.i) 



p j^p I, over J 

Here i labels the different partially conserved axial current (PCAC) mass values for a given Lj. The effect of this 
is to find a y m such that f p for the other values Lj,rriij is as close as possible to the curve obtained from fitting 
Lp,mi^ p . This is summed over all possibilities p. Also, "over" indicates that only i are used such that rriijL^" 1 falls 
within the range of values of rrii iV L v v m , so that the comparison is to an interpolation of the m^pL^" 1 data, rather than 
an extrapolation. Unweighted fits were used so that the approximation to the scaling curve would pass through data 
at small x = mL Vm , where absolute (statistical) errors are largest. (Using a weighted fit reduces our conclusion for 7 
by 4%.) 



Type 




Quadratic 


Co + C\X + C2X 2 


Log quadratic 


Co + ci In x + C2 (In x) 2 


Piece-wise log-linear 


Straight lines connecting data 



TABLE I. Interpolating functions that we use to fit data for 
a fixed L p . In the last case, the straight lines interpolating 
between data are on a semi- log plot. 



For the fitting function we have considered the possi- 
bilities listed in Table HI In the case of the quadratic we 
follow one of the methods of pjl, [23[ . The log quadratic 
fit was motivated by the behavior of the data when L 
is plotted versus \n(mL Vm ), which is close to a parabola. 
The piece- wise log-linear form was used as a third choice 
that trivially passes through the data, giving a reasonable 
interpolation. 



IV. RESULTS 

We have used four observables: the "pion" mass 
77^, the "rho" mass m p , the "ax" mass m ai , and the 
"pion" decay constant f^. These are all obtained from 
standard correlation functions using point sources and 
sinks. We fit the correlation functions with a single 
exponential, allowing the first time ig rst in the fit to 
be large enough for the excited state contributions to 
be negligible. This is determined by looking at the 
mass of the meson as a function of £fi rs t and extract- 
ing the value on the plateau. Five values of bare masses 
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FIG. 1. Collapse of 7r data for y m = 1.46. Here and in the 
other figures, x — mL Vm . 



mod = —1.0, —1.1, —1.165, —1.175, —1.18 on lattices of 
size L/a = 10,12,16,20,24 were simulated, all at /3 = 
2.25. These are the same configurations as were gener- 
ated in [28J , and the values of the PCAC mass and details 
on the simulations are given there. Also note that the size 
of the temporal direction is T — 2L. 

Using these results, and performing the minimization 
described in the previous section, we obtain values for 
y m . In the case of m ai and / w , the quantity Cl/L is 
small, and scaling violations [cf. Eq. (|2.5j) ] can compete 
with the scaling function for small lattices. For this rea- 
son we exclude the small lattices L/a = 10, 12 for these 
channels. The results for y m are summarized in Table 
Hi] It can be seen that each of the channels, and each of 
the fitting methods are consistent with each other within 
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FIG. 2. Collapse of p data for y m — 1.50. 



Observable 


Quadratic 


Log Quad 


PWL 


Combined 


m n 


1.67(93) 


1.26(54) 


1.51(33) 


1.46(27) 


m p 


1.67(88) 


1.37(39) 


1.56(31) 


1.50(23) 


m ai 


1.40(52) 


1.42(27) 


1.41(22) 


1.41(16) 


u 


1.65(22) 


1.49(54) 


1.60(29) 


1.62(17) 



TABLE II. The scaling exponent y m = 1 + 7 for the various 
observables and methods of interpolation. In parentheses, 
jackknife error is shown, obtained from eliminating one rriij 
in all possible ways, in the minimization of H3.1|) . Because we 
use a large number of configurations, O(10 3 ), statistical error 
is negligible by comparison. 



FIG. 3. The 01 data for y m = 1.41. The L = 10, 12 data 
do not fall on the curve, which we interpret as large scaling 
violations. These points were excluded from the fit. 
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errors. The approximate collapse of data is shown in 
Figs. [TJ31 it can be seen in Figs. [3] and 0] that the ex- 
cluded small lattice data does not fall on the scaling 
curve. We interpret this as being due to scaling viola- 
tions, though a thorough study extracting u) would be 
required to demonstrate this. Another interpretation is 
that the theory does not have an IRFP, and so the FSS 
fails for some channels. It is also possible that we are 
seeing the effect of (3 = 2.25 not being close enough to 
the fixed point coupling. However we view the collapse 
seen in Figs. [1] and [2] as favoring our scaling violation 
interpretation. 

We average the twelve values of y m for the four chan- 
nels and three fitting methods, weighted by the jackknife 
errors, to obtain 7 = 0.51. The combined jackknife error 
from the twelve fits is 0.097, and the standard deviation 
of the twelve fits is 0.128. We regard these as two sources 
of systematic error and combine them in quadrature to 
obtain a final estimate of 

7 = 0.51 ±0.16 (4.1) 

In Table IIII1 we compare this estimate with results 
from other groups using a variety of methods. Our re- 
sult is consistent with one of those obtained using the 
Schrodinger functional [llj], perturbation theory, Monte 
Carlo renormalization group and the all-orders hypoth- 
esis of 29]. We obtain a value significantly larger than 



FIG. 4. The data for y m = 1.62. The L = 10, 12 data 
do not fall on the curve, which we interpret as large scaling 
violations. These points were excluded from the fit. 

the one obtained in the FSS studies 0, ^ (1.8a dif- 
ference) and somewhat larger than the Schrodinger func- 
tional study [l7[ (1.3ct difference). 



V. CONCLUSIONS 

We have applied the FSS approach of [23| (one of the 
approaches in [19|]) to MWTC and find values of the crit- 
ical exponent that are in agreement with perturbativc 
results, but somewhat higher than the Schrodinger func- 
tional result with the smallest estimate of error [l7j], by 
1.3cr. While there are significant systematic uncertain- 
ties, which we interpret as being due to scaling violations 
on small volumes, the value of 7 is too small for phe- 
nomenological models of condensate enhancement, which 
requires 7 » 1. The complimentary information obtained 
by the present method suggests that it be applied in other 
gauge theories of interest for conformal or near-conformal 
dynamics. Indeed we expect it to work in any case for 
which the gauge coupling runs very slowly, so that fixed 
point behavior is well approximated on the scales probed 
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Method 


7 


SF [16] 


0.05 < 7 < 0.56 


SF [17] 


0.31 ±0.06 


Perturbative 4-loop [301 


0.500 


Schwinger-Dyson [311 


0.653 


All-orders hypothesis [29] 


0.46 


MCRG [32] 


-0.6 < 7 < 0.6 


FSS [21] 


0.05 < 7 < 0.20 


FSS [22] 


0.22 ± 0.06 


FSS (here) 


0.51 ±0.16 



TABLE III. Summary of all MWTC results for the anomalous 
mass dimension. SF is Schrodinger functional and MCRG is 
Monte Carlo renormalization group. The perturbative result 
7 = 0.500 was also given in the later, corrected version of pjH ]. 
and relied on invariants calculated in [33| . 



by the study that is performed. Unfortunately, as ex- 
plained in [23], a reasonable fit to the FSS assumption 
does not rule in or out the existence of an IRFP, since all 



that is required is a very slow running. 

We have highlighted some of the systematic uncertain- 
ties of the method, and have illustrated how working on 
small volumes hampers the effort to obtain an accurate 
y m . Future work includes simulations on larger volumes 
so that y m can be obtained with greater certainty. Also, 
an improved lattice action should reduce the size of the 
scaling violations, and we are currently working in that 
direction for MWTC and other theories. 



ACKNOWLEDGEMENTS 

This research was supported by the Dept. of Energy, 
Office of Science, Office of High Energy Physics, Grant 
No. DE-FG02-08ER41575. We gratefully acknowledge 
the sustained use of RPI computing resources over the 
course of a year, both the on-campus SUR IBM Blue- 
Gene/L rack, as well as continuous access to 1-4 racks 
of the sixteen IBM BlueGene/L's situated at the Com- 
putational Center for Nanotechnology Innovation. We 
thank Tom DeGrand for extensive discussions and com- 
ments. We also thank Luigi Del Debbio, Julius Kuti and 
Francesco Sannino for helpful comments. 



[1] L. Susskind, Phys. Rev. D20, 2619 (1979) 

[2] S. Weinber g, |Phys. Rev. D19, 1277 (197971 

[3] B. Holdom, [Phys. Rev . D24, 1441 (1981) 

[4] B. Holdom, Phys . Lett. B150, 301 (1985) 

[5] K. Yamawa ki, M. Bando, and K.-i 

|Phys. Rev. Lett. 56, 1335 (1986)| 

[6] M. Bando, K. Matumoto , 

[7] ~ 



Matumoto, 



and K. Yamawaki, 



s. Lett. B178, 308 (1986)] 
T. W. Appelquist, D. Karabali, and L. C. R. Wijeward- 
hana, |Phys. Rev. Lett. 57, 957 (1986)| 
T. Appelquist and L. C. R. Wijewardhana, 



|Phys. Rev. D35, 774 (1987)] 
T. Appelquist and LT C. 



R. Wijewardhana, 



[8] 
[9] 

|Phys. Rev. D36, 568 (19~87T] 

[10] E. Eichten and K. D. Lane, |Phys. Lett. B90, 125 (1980)| 
[11] K. D. Lane and E. Eichten, 

|Phys. Lett. B222, 274 (1989) 
[12] F. 



Sannino 



and 



K. 



Tuominen, 



Phys. Rev. D71, 051901 (2005), arXiv:hep-ph/0405209 



[13] L. Del Debbio, PoS LATTICE2010, 004 (2010). 
[14] H. S. 



Fukano 



and 



San- 



[15] 

[16] 

[17] 
[18] 



arXiv: 1005.3340 [hep-ph 



Phys.Rev. D82, 035021 (2010) | 



~Y. Shamir, K Svetitsky, and T. De 

Grand, | Phys. Rev. D78, 031502 (2008) 



arXiv:0803.1707 [hep-lat 



F. Bursa, L. Del Debbio, L. Keegan, C. Pica, 
and T. Pickup, iPhys.Rev. D81, 014505 (2010) | 
|arXiv:0910.4535 [hep-ph]j ~ 

T. DeGrand, Y. Shamir , and B. Svetitsky, (2011), 
|arXiv:1102.2843 [hep-lat]] 

Z. Fodor, K. Holland, J. Kuti, D. Nogradi, and 
C. Schroeder, PoS LATTICE2008, 058 (2008), 



[19 
[20 

[21 

[22 

[23 
[24 

[25 

[26 
[27 
[28 

[29 
[30 
[31 



|arXiv:0809.4888 [hep-lat] 




T. DeGrand, J 


3 hys.Rev. D80, 114507 (2009) | 


|arXiv:0910.3072 [hep-lat] 




L. Del Debbio and R. Zwicky, 



Phys.Rev. D82, 014502 (2010) , 
|arXiv: 1005.2371 [hep-phjj ~ 

L. Del Debbio, B. Lucini, A. Patella, C. Pica, 
and A. Rago, |Phys.Rev. D82, 014509 (2010) 
|arXiv: 1004.3197 [hep-lat 



and A. Rago, |Phys.Rev. D82, 014510 (2010) 


|arXiv: 1004.3206 [hep-lat] 




T. DeGrand, 


Phys.Rev. D84, 116901 (2011) 


|arXiv:1109.1237 [hep-lat] 





Z. Fodor, K. Holland 
C. Schroeder, et al, 
|arXiv:1104.3124 [hep-lat 
T. Appelquist^ G~ Fleming, M. Lin 
and 



J. Kuti, D. Nogradi, 
Phys.Lett. B703, 348 (2011)1 



E. Neil, 



D. Schaich, Phys.Rev. D84, 054501 (2011) 



|arXiv: 1106.2148 [hep-lat]p 
J. Kuti, Talk given at Lattice 2011, July 10-16, 2011, 
Squaw Valley, USA. 
S. M. Bhattacharjee and 
| Jour, of Phys. A34, 6375 (2001) | 
J. Giedt and 



Seno, 



E. 



Weinberg, 



Phys.Rev. D84, 074501 (2011) 
|arXiv:1105.0607 [hep-lat 1 1 

C. Pica and F. Sannino, Phys.Rev. D83, 116001 (2011) 



|arXiv: 1011.3832 [he p-ph 

C. Pica and F. Sannino, IPhys.Rev. D83, 035013 (2011) 
|arXiv: 1011.5917 [hep-ph]] - 

T. A. Ryttov and R. Shrock, 

IPhys.Rev. D83, 056011 (2011)| 



|arXiv:1011.4542 [hep-ph]| 
[32] S. Catterall, L. Del Debbio, J. Giedt, and L. Keegan, 

(2011), |arXiv:1108 3794 [hep-ph]) 
[33] M. Mojaza, C. Pica, and F. San- 

nino, |Phys.Rev. D82, 116009 (2010) | 



arXiv: 1010.4798 [hep-ph]] 



